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ABSTRACT 

Intermediate resolution spectroscopy from the ESO Very Large Telescope is analysed 
for 63 photometrically selected low-mass (0.08-0.30 M ) candidates of the open cluster 
NGC 2547. We have confirmed membership for most of these stars using radial veloci- 
ties, and found that lithium remains undepleted for cluster stars with / > 17.54 ±0.14 
and K s > 14.86 ± 0.12. From these results, several pre-main-sequence evolutionary 
models give almost model independent ages of 34-36 Myr, with a precision of 10 per 
cent. These ages are only slightly larger than the ages of 25-35 (±5) Myr obtained using 
the same models to fit isochrones to higher mass stars descending towards the zero 
age main sequence (ZAMS), both in empirically calibrated and theoretical colour- 
magnitude diagrams. This agreement between age determinations in different mass 
ranges is an excellent test of the current generation of low-mass pre-main sequence 
stellar models and lends confidence to ages determined with either method between 
30 and 120 Myr. 

Key words: stars: abundances - stars: late-type - open clusters and associations: 
individual: NGC 2547 



1 INTRODUCTION 



Absolute ages in young open clusters and star forming re- 
gions are difficult to determine and usually very model 
dependent, yet are crucial to solving many stellar evolu- 
tion problems - from the lifetimes of protoplanetary disks 
and progression of angular momentum loss, to calibrating 
cooling models and the initial-final mass relation for white 
dwarfs. Nuclear turn-off ages from high mass stars may be 
uncertain by factors of two, depending on the treatment of 
core convection, rotation and mass loss (Chiosi, Bertelli & 
Bressan 1992; Meynet et al. 1993; Meynet & Maeder 1997, 
2000). Ages from fitting isochrones to cool, low-mass stars, 
descending towards the ZAMS may be equally inaccurate: 
they depend on detailed modelling of stellar atmospheres, in- 
cluding molecular opacities, convection and turbulence (see 
Baraffe et al. 1998, 2002) and observational comparisons suf- 
fer from uncertainties in the conversion from model parame- 
ters (Lboi and T c ff) to magnitudes and colours (e.g. Stauffcr, 
Hartmann & Barrado y Navascues 1995). 

The Lithium Depletion Boundary (LDB) is a new, im- 



* Based on observations collected with the VLT/UT2 Kueyen 
telescope (Paranal Observatory, ESO, Chile) using the 
FLAMES/GIRAFFE spectrograph (Observing run 072.D- 
0406 A). 



portant, and possibly less model dependent method for de- 
termining the ages of young clusters. The initial Li content 
of a low-mass star begins to burn in p, a reactions, as stars 
contract along PMS tracks and core temperatures reach 
2.5 x 10 6 K. Below O.O65M0, Li is never burned because 
the core does not become hot enough. At higher masses, 
the time taken to reach Li-burning temperatures is a sen- 
sitive function of mass and hence a very sensitive function 
of luminosity (Bildsten et al. 1997; Ushomirsky et al. 1998). 
Low-mass stars approaching the ZAMS are fully convective, 
so their entire Li content is then rapidly exhausted (in less 
than a few Myr after Li burning is initiated) and thus in 
a cluster, the luminosity at which photospheric Li remains 
undepleted can give a very precise age estimate and is most 
sensitive for ages of 10-200 Myr. 

The small experimental uncertainties in the LDB 
method are due to the conversion from optical and near in- 
frared magnitudes to luminosities via cluster distances and 
empirical bolometric corrections, and the difficulty in locat- 
ing the LDB among a population of very faint, low-mass 
stars (see Jeffries & Naylor 2001). Systematic errors are 
likely to be quite small. Remarkably, different evolutionary 
models and even analytic calculations yield almost identical 
(to ±10%) ages, despite very different treatments of atmo- 
spheres, convection, opacities and equations of state (see for 
example Burke, Pinsonneault & Sills 2004). It is this model 
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independence which makes it vital to measure the LDB in as 
many young clusters as possible. Such clusters would define 
empirical isochrones that could be used to link together and 
calibrate the uncertain physics in both high- and low-mass 
stellar evolution models. 

LDB ages have so far only been determined for the 
Pleiades (125±8 Myr), Alpha Per (90±10Myr) and IC2391 
(53 ± 5 Myr) clusters (Stauffer, Schultz & Kirkpatrick 1998; 
Stauffer et al. 1999; Barrado y Navascues, Stauffer & Pat- 
ten 1999; Barrado y Navascues, Stauffer & Jayawardhana 
2004). These ages are older by factors of ~ 1.5 than nu- 
clear turn-off ages determined from high mass evolutionary 
models with zero convective overshoot. The fractional dis- 
crepancy may become smaller towards younger ages (and 
hence higher masses at the turn-off), hinting that additional 
convective overshoot is mass-dependent. 

For IC2391, the LDB age is also about 15 Myr older 
than the age indicated by fitting isochrones to the low-mass 
(0.3 < M < 1.2M©) PMS cluster members. This may also 
be true for the cluster NGC 2547, which has a nuclear turn- 
off age of 55 ± 25 Myr (Jeffries & Tolley 1998), a low-mass 
isochronal age of 20-35 Myr (Naylor et al. 2002) a distance 
modulus of 8.10±0.10 (Naylor et al. 2002) and a reddening, 
E(B -V) = 0.06 ± 0.02 (Claria 1982). A previous attempt 
to detect the LDB in NGC 2547 was made by Oliveira et 
al. (2003), but resulted in lower limits to the cluster age. No 
Li-rich objects were found for / < 17.2 and by considering 
the average spectrum of several cluster candidates there was 
some evidence that the LDB lay between 17.8 < I < 18.3. 
This suggested that the LDB age of the cluster was certainly 
more than 30 Myr and probably more like 38-46 Myr and 
hence older than the low-mass isochronal age. Oliveira et al. 
suggested that this discrepancy might indicate shortcomings 
in the commonly used low-mass evolutionary models. 

In this paper we present optical spectroscopy from the 
ESO Very Large Telescope (VLT) that unambiguously iden- 
tifies a substantial population (~ 50) of very low-mass stars 
in NGC 2547 and now precisely determines the location of 
the LDB. In section 2 we discuss the selection of targets, 
the fibre spectroscopy performed at the VLT and the data 
reduction process. In section 3 we report measurements of 
relative radial velocities and the Li I 6708A feature in the 
cluster candidates, and our assessment of cluster member- 
ship. Section 4 deals with the location of the LDB and the 
age of NGC 2547. The results are discussed in section 5 and 
conclusions drawn in section 6. 



2 VLT FIBRE SPECTROSCOPY 
2.1 Target selection 

The targets for VLT fibre spectroscopy were chosen from 
the RIZ survey of Jeffries et al. (2004) . This catalogue con- 
tains precise astrometric and photometric information for 
stellar objects in 0.855 square degrees around NGC 2547 
and is substantially complete to / = 19.5. Fifty four tar- 
gets were chosen from the lists of candidate cluster members 
with 16.4 < I < 19.1 which were in turn selected from their 
photometric colours and proximity to empirically calibrated 
30 Myr isochrones generated from the models of D'Antona & 
Mazzitelli (1997) and Baraffe et al. (2002) (see Jeffries et al. 



2004 for details). Using an intrinsic distance modulus of 8.1 
and extinction, Ai = 0.112 (Rieke & Lebofsky 1985), clus- 
ter members in this magnitude range will have approximate 
masses of 0.08-0.30 Mq according to the models of Baraffe 
et al. (2002). 

These targets are only a small, representative subset of a 
larger population of NGC 2547 cluster members in this mag- 
nitude range, that is spread over more than a square degree. 
A further nine objects were observed that were slightly more 
distant from the cluster isochrones and not classed as mem- 
bers by Jeffries et al., but which were also contained within 
the 25 arcminute diameter VLT field of view. The precision 
of the photometry is about ±0.04 in / and ±0.07 in R — I 
for the faintest object in the sample (with / = 19.05), but 
is more like ±(0.01 — 0.02) in colour and magnitude for the 
majority. 



2.2 Observations 

The targets were observed on 19 and 20 January 2004 us- 
ing the FLAMES instrument mounted on the VLT-Kueyen 
(UT2) 8.2-m telescope. The GIRAFFE spectrograph and 
MEDUSA fibre system were used in combination with the 
600 lines mm -1 grating and an order sorting filter to give 
spectra with a resolving power of about 8000 over the wave- 
length range 6438A < A < 7184A. The MEDUSA system 
has 130 fibres, each with a projected diameter of 1.2 arc- 
seconds on the sky, that can be positioned within a circular 
area of diameter 25 arcminutes. 

The targets were observed in 8 separate observing 
blocks (OBs), each of which yielded 2589 s of observing time 
on sky and all of which were centred at RA = 08h 10m 06.6s, 
Dec= — 49dl5m43s. A number of fibre configurations were 
used, such that most objects (and all with / > 18.5) were 
observed in all 8 OBs, a few targets were observed in 5 or 6 
OBs and one target (number 35 - see Table0 was observed 
in only 2 OBs. This was done to maximise the number of 
targets observed at the required level of precision. Six of the 
OBs were obtained on 19 January 2004 between 03:14 UT 
and 08:43 UT. The other two OBs were obtained on 20 Jan- 
uary 2004 between 03:09 UT and 07:28 UT. About 50 fibres 
which could not be allocated to useful targets were placed on 
blank sky regions in each configuration. The seeing measured 
at the telescope was 0.5-0.9 arcseconds and the moon was 
below the horizon during the observations. Flat-field, bias 
and Thorium- Argon arc lamp calibrations were performed 
at the start of each observing night. 



2.3 Data Reduction 

The calibrations and observational data were reduced using 
version 1.10 of the GIRAFFE pipeline software (Blecha et 
al. 2003). Master bias and flat-field calibrations were pro- 
duced for each night along with fits tables containing the 
locations of the fibre spectra (from the master flat-field) on 
the 2Kx4K EEV CCD detector and wavelength calibration 
solutions for all the fibres taken from the arc lamp spec- 
tra. Each target exposure was bias-subtracted, and scattered 
light was removed by modelling the remaining signal be- 
tween the individual spectra. Target and blank sky spectra 
were extracted using an optimal algorithm. Pixel-to-pixel 
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and fibre-to-fibre sensitivity variations were corrected for by 
dividing by a normalised set of spectra extracted from the 
master fiat-field using the same apertures. Note that no at- 
tempt is made to flux calibrate the spectra in an absolute 
sense, but that most of the instrumental response is removed 
by dividing by the flat-field continuum. The corrected spec- 
tra were then rebinned to a wavelength scale of 0.2A per 
pixel, slightly larger than the average wavelength coverage 
of a single detector pixel. 

Sky subtraction was achieved by taking the median of 
groups of five blank sky spectra and averaging these me- 
dians to form a master sky spectrum which could be sub- 
tracted from each target spectrum. The standard deviation 
among the continuum levels in the individual sky spectra 
was about 7 per cent (1-sigma). Part of this arises from 
uncertainties in the fibre-to-fibre sensitivity corrections, re- 
sulting largely from non-uniform illumination of the fibres 
by the flat-field lamp (L. Pasquini, ESO, private commu- 
nication). There may also be a contribution from scattered 
continuum radiation from a bright, spatially varying, Hn re- 
gion which coincides with the NGC 2547 line of sight and 
is clearly visible in 7?-band images. The Hn region is even 
more evident in the bright night-sky emission lines, where 
variations of 20-30 per cent are found even between sky fi- 
bres situated only a few arcminutes apart. 

Uncertainties in the subtracted sky continuum are of 
no concern for this paper as far as measuring the continuum 
level close to the LiI6708A feature. The sky continuum is 
just equal to the signal from the targets at 7 ~ 19 and 
hence for an unresolved line of width ~ 0.8A there may be 
an additional equivalent width error of about 0.05A, which 
is smaller than the statistical errors in the equivalent width 
measurements even for the faintest targets (see section 3.3). 
However, the strong sky lines are up to 200 times stronger 
than the sky continuum and so sky-subtraction in the imme- 
diate vicinity of these wavelengths (Ha, Sn 6717A, 673lA 
etc.) was more problematic. 

Multiple spectra for targets were combined using an it- 
erative technique that rejected pixels deviating by more than 
5-sigma from from the average of pixels at the same wave- 
length in all the normalised spectra of that target. Three 
separate combinations were produced: all of the spectra; all 
spectra from 19 January 2004; and all spectra from 20 Jan- 
uary 2004. Where only two spectra of an object existed in 
the combination, these spectra were simply summed. The 
theoretical levels of signal-to-noise ratio (SNR) achieved in 
the averaged spectra from all the target exposures, based on 
measurements of the CCD gain and readout noise and prop- 
agating errors in the bias subtraction, flat-fielding, scattered 
light subtraction and sky subtraction, ranged from 6 to more 
than 50 per 0.2A pixel in the faintest and brightest targets 
respectively. A sample of our spectra covering the full range 
of wavelengths, colours and magnitudes is shown in Fig. 



3 ANALYSIS 

3.1 Spectral Indices 

To check that all the observed targets had spectral types 
appropriate for their colours, we calculated the TiO (7140A) 
and CaH (6975A) narrow band spectral indices (see Oliveira 



et al. 2003) . These indices are temperature sensitive and we 
expect a smooth relationship between their values and the 
R — I colour index of NGC 2547 cluster members. 

None of our targets are found to deviate significantly 
from the trend defined by the bulk of the data. How- 
ever, this is not unexpected. Any contamination in our 
sample is certain to be dominated by foreground field M- 
dwarfs with almost identical narrow band indices at the 
same colour (Jeffries et al. 2004). The CaH index does have 
some sensitivity to gravity but this is unlikely to distinguish 
between cool pre-main-sequence objects at 30-50 Myr and 
main-sequence dwarfs. Indeed, targets which are deemed 
cluster non-members by virtue of their relative radial ve- 
locities (see section 13.21 are indistinguishable from cluster 
members on the basis of these narrow band indices. Table 1 
lists spectral types for our targets, estimated by comparing 
their TiO (7140A) narrow band indices with the those mea- 
sured for stars with known spectral type (see Oliveira et al. 
2003). We estimate these spectral types are good to ± half 
a subclass. 

3.2 Radial velocities 

Although we were unable to obtain observations of radial 
velocity standard stars during our observations, the relative 
radial velocities (RVs) of our targets proved a useful tool to 
discriminate cluster members from non-members and also 
identify short period (of order 10 days or less) binary systems 
in those targets which show other indications of membership 
(i.e. a detection of Li - see section 13.31 . Members of the 
cluster (bar short period binaries) should have RVs with a 
dispersion of order 1 km s _1 (Jeffries, Totten & James 2000), 
so given that we expect the bulk of our targets to be cluster 
members, selection in a narrow range of relative RVs should 
exclude the vast majority of non-members. 

Radial velocities were calculated by cross-correlating 
the region 6950 < A < 7150A which includes prominent 
molecular band heads and only a few very weak nebular 
emission lines and telluric features. We chose one of the tar- 
gets (number 30 - see Table with a good SNR and an 
intermediate colour and spectral type to act as the RV tem- 
plate for all the other objects. All the spectra from each 
night were cross-correlated with the spectrum of object 30 
obtained on 19 January 2004. Hence most objects have two 
relative RV measurements. 

The average relative RV for the cluster members and 
the precision of the RV measurements were estimated by 
considering the subsample of targets that were Li-rich and 
hence almost certain cluster members (see section 13.51 . 
but excluding three objects which had clearly discrepant 
or varying RVs from night to night. We found an aver- 
age RV of (—0.9 ± 0.3) kms -1 with a standard deviation 
of 1.7 kms -1 on 19 January 2004 and an average RV of 
(+0.0 ± 0.6) kms _1 with a standard deviation of 2.6kms~ 1 
on 20 January 2004. As the combined data from 19 Jan- 
uary had a higher SNR, the smaller standard deviation is 
expected. 

Cluster members were required to have RVs within na 
of the average RV of the Li-rich objects on 19 January 2004 
and a consistent RV measured on 20 January 2004. The 
value of n was adjusted depending on the empirical SNR 
(defined in the next section) of the spectrum in question. For 
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Figure 1. A representative sample of our reduced spectra covering the full range of colours and magnitudes. The spectra are normalised 
to unity at 6705 A and offset. The labels refer to the identifiers in Tabic HI 



SNR > 20, we used n = 2, for 10 < SNR < 20, n = 3 and for 
SNR < 10, n = 4. Before comparing the RVs from 19 and 
20 January 2004, the mean discrepancy of 0.9 km s -1 was 
subtracted from the data taken on 20 January. The relative 
RVs for the targets on each night are listed in Table Q 

This process yielded 52 objects deemed to be cluster 
members, 9 objects deemed to be RV non-members and a 
further 2 objects where the RV appeared to be variable be- 
tween data from the two nights. For some targets the cross- 
correlation function was markedly broadened, perhaps as a 
result of rapid rotation (see below) or binarity. These RV 
measurements are indicated in Tabled 



3.3 Lithium 

The main aim of our experiment was to search for Li in 
the form of the LH6708A resonance doublet. This feature 
should be strong in cool stars with undepleted Li - with 
an equivalent width (EW) of 0.5A to 0.6A according to the 



curves of growth presented by Zapatero-Osorio et al. (2002). 
Figure |21 shows a selection of our spectra in this wavelength 
region, where data combined from both nights are shown. 

The EW of the 6708A feature was estimated with re- 
spect to a pseudo-continuum defined by a linear fit to the 
data in the wavelength ranges 6700 to 6705A and 6710 to 
6713A. The upper limit to this latter continuum region was 
chosen to coincide with a downturn in the spectra due to 
a molecular bandhead, but also to avoid regions within a 
few angstroms of the Sn 6717A sky line, which was poorly 
subtracted in many of our spectra and contained very noisy 
data in any case. The residuals to this linear fit gave an 
empirical estimate of the SNR which was used in turn to 
estimate either an uncertainty in the Li EW or a 2-sigma 
upper limit where no line was seen. The empirical SNR is 
bound to underestimate the true SNR of the spectra because 
there are molecular features in the pseudo-continuum that 
become stronger at cooler temperatures. Hence our EW un- 
certainty estimates and upper limits should be conservative 



© 2004 RAS, MNRAS 000.HHT71 



The Lithium depletion boundary of NGC 2547 5 



6.0 - 



4.0 - 



3.0 - 



2.0 - 



1.0 



1 


1 1 


5 


/ 














27 

i 


l l 



- 6.0 - 



- 5.0 - 



- 4.0 - 



- 3.0 - 



- 2.0 - 



1.0 



1 


i i 












A ^ 






37 kkt a 








i 


1 1 



- 6.0 - 



- 5.0 - 




6690 6700 6710 6720 
Wavelength (A) 



67306690 6700 6710 6720 
Wavelength (A) 



67306690 6700 6710 6720 
Wavelength (A) 



6730 



Figure 2. Representative spectra in the vicinity of the Li I 6708A spectral feature. The spectra are labelled according to the identifiers 
in Table H\ and plotted in 0.4A bins. The gaps in the spectra are the location of the prominent S II line found in the sky spectra, which 
could not be accurately subtracted. The spectra are normalised to unity around 6705A and offset. The vertical dashed line indicates the 
expected position of the Li I 6708A feature for members of NGC 2547. 



and our EWs for the coolest stars may be systematically 
underestimated. 

The EWs were estimated by fitting and integrating a 
Gaussian absorption feature beneath the linear continuum. 
Uncertainties were estimated according to the approximate 
formula 



AEW = 



VFWHM x p 
SNR 



(1) 



where FWHM and p are the FWHM of the Gaussian and the 
size of the pixels in A. The Li EWs or 2-sigma upper limits 
are listed in TablcQ]along with the empirical SNR estimate. 
For target 63, which showed a large RV shift between the 
two nights, the quoted results are from the first night. 

Most of our targets have FWHM in the range 0.7 to 
TOA as expected from the resolving power of the spec- 
trograph. A few objects showed Li lines which were much 
broader than this, with FWHM of 1.5-2.2A. This is likely 



due to rapid rotation either in single objects simply because 
they are young (a phenomenon certainly seen in higher mass 
objects in the cluster - Jeffries et al. 2000) or because they 
are part of close, tidally locked binary systems. These broad- 
ened objects are indicated in TableQ The data are of insuffi- 
cient quality to estimate projected rotational velocities, but 
they must be of order 50kms _1 to produce the broadening 
observed. 
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Table 1. Targets observed in this paper. The columns list identifiers used in this paper, identifiers from the catalogue of Jeffries et al. 
(2004) a "N" following the name indicates it was not considered a photometric member in that paper, positions, I magnitudes, R — I 
and I — Z colours from Jeffries et al. (2004), JHK S magnitudes from 2MASS (Cutri et al. 2003), estimated spectral type (see section 
3.1), relative radial velocity (nights 1 and 2, a "b" indicates a broadened cross-correlation function — see section [3.2l . whether the radial 
velocity appears variable (V), the empirical SNR of the spectra around the Li 6708A feature, the measured equivalent widths for this 
feature or a 2-sigma upper limit (a "b" indicates the line is broadened — see section [3.31 and our membership assessment status for the 
object - (1) radial velocity member with no Li, (2) radial velocity member with Li, (3) radial velocity non-member with Li, (4) radial 
velocity non-member with no Li. Photometric uncertainties are 0.01 to 0.04 mag in / and I — Z, 0.01 to 0.06 mag in R — I and 0.04 to 
0.20 mag in JHK S (from brightest to faintest). 
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3.4 Ha emission 

Ha emission was seen in almost all our sky-subtracted spec- 
tra and in principle could be used as an additional member- 
ship criterion (e.g. Barrado y Navascues et al. 2004). How- 
ever, whilst chromospheric Ha emission is expected from 
NGC 2547 members it would also be expected from a large 
fraction of any possible field M-dwarf contaminants with 
similar spectral types (Gizis, Reid & Hawley 2002). Unfor- 
tunately the presence and strength of Ha emission could 
not be accurately measured because the Ha emission line 
present in the sky spectrum is much stronger. The peak in- 
tensity of the Ha in the mean sky spectrum is 20-200 times 
the continuum intensity around Ha in our targets. Further- 
more, the strength of the sky Ha emission varies by 20-30 
per cent from fibre to fibre because of the background H n 
region (see section 2.3), resulting in effective EW errors of a 
few Angstroms even in the brightest stars of our sample. Slit 
spectroscopy will be required to obtain accurate estimates 
of chromospheric Ha emission in these stars. 

What we were able to do was search for extended Ha 
emission beyond the wings of the sky line. In particular, 
White & Basri (2003) have shown that emission lines with 
a full width at 10 per cent of maximum of more than 
270 km s -1 are indicative of ongoing accretion activity. Such 
objects would also typically have emission EWs of more than 
20A and hence extended emission wings easily visible in 
our spectra. Following the method detailed in Kenyon et 
al. (2004) we simulated a Gaussian Ha emission line with 
an EW of 20A and full width at 10 per cent of maximum of 
270 kms -1 for each of our spectra. No emission beyond this 
simulated profile (appropriately broadened for those objects 
thought to be rapid rotators - see section 3.3) was found. 
This result is not surprising. Jeffries et al. (2000) found no 
accretion signatures in higher mass objects in NGC 2547 and 
a recent Spitzer survey of the cluster by Young et al. (2004) 
finds less than 7 per cent of the low-mass cluster members 
have inner dust discs betrayed by an infrared excess. Nev- 
ertheless, combined with the precise LDB age derived for 
the cluster in this paper, NGC 2547 provides a useful upper 
limit to the lifetime of this early phase of stellar evolution. 



3.5 Cluster membership 

The measurements discussed above allow an assessment of 
cluster membership. First, we assume that the 23 objects 
showing an Li feature with an EW > 0.2 A and an RV con- 
sistent with cluster membership are cluster members. The 
rationale is that curves of growth presented by Zapatero- 
Osorio et al. (2002) suggest that the Li abundance is de- 
pleted by about a factor of 100 or less in these objects (see 
section 0J, which is also the criterion we will adopt for judg- 
ing the location of the LDB. Field M dwarfs are extremely 
unlikely to show an Li feature of this strength. 

In addition there are two stars which show Li at this 
level but have an RV inconsistent with cluster membership. 
One of these (star 60) is probably a short period binary 
member of the cluster, the other (star 5) is unlikely to be 
a cluster member based on its colour and magnitude (see 
section 16.21 . The 9 objects which show no Li and which 
have an RV inconsistent with cluster membership could be 
cluster binaries though this is unlikely. Five of these objects 




1.5 1.6 1.7 1.8 1.9 2 2.1 2.2 

R-I 

Figure 3. Li EWs versus R — I colour for Li-rich and Li-poor 
RV members of NGC 2547. Downward pointing triangles indicate 
2-sigma upper limits. Encircled points are those with significant 
rotational broadening — see section 3.3. 

are fainter than our estimate of the LDB (see sectionQjJ and 
should have exhibited Li if they were members and some 
of the others look unlikely to be cluster members on the 
basis of their colours and magnitudes (see sect ion lo"21 . Until 
there is further evidence we classify all these objects as non- 
members. 

Finally, there are 29 objects with an RV consistent with 
cluster membership but showing no Li. These are very likely 
to be cluster members, although given the ~ 10 km s _1 range 
over which cluster members are selected, it is possible that 
one or two contaminating non-members remain in this sub- 
set. There are unlikely to be more than this because we see 
no Li-poor, RV members fainter than the location of the 
LDB. 



4 LITHIUM ABUNDANCES 

We have estimated Li abundances, .A (Li) (where -A(Li) = 
log[A r (Li)/7V(H)] + 12), for those objects which are RV mem- 
bers of the cluster. We use the curves of growth for cool 
stars presented by Zapatero-Osorio et al. (2002), and ex- 
tended to A(Li)= 0.0 by Pavelenko & Zapatero-Osorio (pri- 
vate communication) , in order to convert our EW values into 
Li abundances. The difficulties in doing this are described 
extensively in Jeffries et al. (2003), the main problem being 
that the definition of the pseudo-continuum used to derive 
the EW is not necessarily consistent between objects mea- 
sured at different spectral resolutions or with different de- 
grees of rotational broadening. It is also unknown to what 
extent there are model-dependent systematic errors in the 
cool atmospheres used to determine the curves of growth. 
For these reasons the absolute values of the Li abundances 
should be treated with caution. 

To estimate T e s values for our objects we have used 
an empirical relationship between intrinsic R — I and T e ff 
defined by making a 120 Myr isochrone from the Chabrier 
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Figure 4. Li abundances versus T c g for RV members of NGC 
2547 from this paper (open symbols, encircled points indicate 
those with significant rotational broadening — see section 3.3) and 
from Jeffries et al. (2003 - filled symbols). Downward pointing 
triangles indicate 2-sigma upper limits. For comparison we show 
model predictions for the progression of Li depletion at 30 and 
50Myr, taken from the models of Chabrier & Baraffe (1997). 



& Baraffe (1997) evolutionary models fit published Pleiades 
RI data at a distance of 132 pc (see section 6.3.2 for details). 
This relationship is shown in Fig. 1101 and the T e g values for 
the cluster members (derived assuming that E(R — I) = 
0.043) are given in Table [5] Statistical uncertainties which 
arise from uncertainties in the R — I photometry are small. 
However to allow for some variability we assign uncertainties 
of ±100 K for the purpose of estimating possible Li abun- 
dance errors. The Li abundances were then obtained using 
a bicubic spline interpolation of the curves of growth. Many 
of our EW upper limits imply Li abundances that fall well 
below the lowest value in our grid, but we conservatively 
assign an upper limit of A(Li)< 0.0 to these. The statistical 
errors in the abundances were found by perturbing the EW 
and T e ff by their error bars and combining the two sources 
of error in quadrature. The EW uncertainty is dominant be- 
cause the curves of growth become nearly independent of 
Teff for Teff < 3500 K. For this reason, systematic uncertain- 
ties in the T e ff scale, which are quite likely to be of order 
100-200 K at these cool temperatures, are not important. Li 
abundances and uncertainties are listed in Tabled 

In Figs |21 and [I] we show plots of the Li EWs versus 
R — I colour and Li abundances versus T e g for the NGC 
2547 RV members. In the latter plot we have also included 
Li abundances for higher mass objects in the cluster derived 
using the same curves of growth (from Jeffries et al. 2003). 
A number of points arise from these two figures. 

(1) There is a clear pattern of undetectable Li for 
warmer stars (R — I < 1.75). At slightly cooler tempera- 
tures there exist stars both with and without detectable Li 
and then for R — I > 1.84 only Li-rich objects are found. 

(2) The overall pattern of Li depletion agrees qualita- 
tively (but not quantitatively) with those predicted by PMS 
evolutionary models (see also Fig. 9 and the discussion in 



Jeffries et al. 2003). In particular the sharp upturn at cooler 
temperatures - the LDB - is present. However it is impos- 
sible to define a precise or accurate age for NGC 2547 on 
the basis of Fig. 21 because: (i) objects with and without 
Li co-exist at the crucial temperatures because of scatter 
caused by photometric errors, variability and unresolved bi- 
narity; (ii) the T c g at the LDB is model dependent by of 
order ±100 K (see Table |3J, which translates to significant 
age uncertainties; (iii) there are systematic uncertainties of 
a similar size in the conversion from R — I to T c g. These 
problems are solved in the next section by considering the 
data in the Hertzsprung-Russell diagram and finding the lu- 
minosity of the LDB. 

(3) There appears to be a scatter in the Li abundances 
at low temperatures and even among the subset of objects 
with luminosities well below the LDB. This scatter arises 
directly from a spread in the EWs at a given colour. We 
believe that the spread in the EWs is real. Although some 
readers may question whether the approximation of equa- 
tion 1 is strictly valid, it is unlikely to underestimate the 
EW errors by the factor of three that would be required 
to make them consistent with no spread. Systematic EW 
uncertainties are likely to be a monotonic function of T e ff 
and should not introduce scatter at a given colour. A simi- 
lar Li EW scatter has been observed by Zapatero-Osorio et 
al. (2002) and Kenyon et al. (2004) for very low mass stars 
with T ff < 3500 K in the young a Orionis cluster and by 
Barrado y Navascues et al. (2004) in IC2391, but the spec- 
tra here are of higher quality and more importantly, taken 
at a homogeneous and relatively high resolution. 

Whether the spread in abundances is real is a different 
matter. There is no theoretical explanation of such a scat- 
ter. All evolutionary models predict a sharp change from no 
Li to essentially cosmic Li over a small temperature range. 
There should be few, if any, transitional objects 1 . Instead 
we think this problem may be related to our understanding 
of the upper atmosphere of very cool stars. Because the Li 
feature is saturated in these cool dwarfs it is very suscep- 
tible to NLTE effects or perhaps overionization caused by 
an overlying chromosphere (e.g. Houdebine & Doyle 1995). 
This would cause a weakening of the Li feature in the most 
active stars and it is notable that: (i) the peak levels of 
Li abundance are roughly that expected for the cosmic Li 
abundance A(Li)~ 3.3; (ii) the rapid rotators tend to have 
lower Li EWs and abundances, although it would be difficult 
to discount problems with blending and continuum defini- 
tion as an explanation for this. Given slit spectroscopy (to 
overcome the Ha sky subtraction problems) or very deep X- 
ray observations it should be possible to look for an inverse 
correlation between Li strength and chromospheric/coronal 
activity. 

(4) Despite the scatter, the objects do still clearly bi- 
furcate between those objects with Li (with an EW> 0.2A) 
and those with much lower upper limits. This division cor- 
responds roughly to where 99 per cent of the original Li has 



1 Target number 21 with a weak but detectable Li line and 
A(Li)= 0.0 ± 0.3 may be an example of a transitional object, 
but could also be an example of an unresolved binary system 
with one component above and one component below the LDB. 
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Table 2. Effective temperatures derived from the intrinsic R — 
I colour and Li abundances (where A(Li)= log[A r (Li)/A r (H)] + 
12) obtained from the curves of growth of Zapatero-Osorio et al. 
(2002) for RV members of NGC 2547 (status 1 and 2 objects from 
Table Columns 1 and 4 are the identifiers used in Table H\ 



ID 



(K) 



A(Li) 



ID 



(K) 



1 


3272 


< 0.0 


31 


3149 


2 


3328 


< 0.0 


32 


3184 


1 


3423 


< 0.0 


34 


3374 


6 


3361 


< 0.0 


35 


3246 


7 


3291 


< 0.0 


36 


3160 


8 


3354 


< 0.0 


37 


3299 


9 


3386 


< 0.0 


38 


3322 


10 


3442 


< 0.0 


39 


3189 


11 


3348 


< 0.0 


41 


3099 


12 


3127 


< 0.0 


42 


3214 


13 


3402 


< 0.0 


11 


3251 


11 


3282 


< 0.0 


45 


3176 
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19 


3003 


19 
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20 
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21 
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< 0.0 


54 
55 


3217 
3040 


25 


3275 
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Figure 5. I,R-I colour-magnitude diagram for our targets. Open 
squares represent objects with an RV consistent with cluster 
membership but which have no Li, filled squares are those RV 
members with Li, small filled squares are Li-rich stars with an RV 
discrepant from the cluster and crosses represent objects with dis- 
crepant velocities and no Li. Encircled points represent those ob- 
jects with discernible rotational broadening (see section 3.3). The 
solid line is an empirically calibrated (see text) 30 Myr isochrone 
generated from the models of Baraffe et al. (2002), the dashed line 
is a similar isochrone generated from the D'Antona & Mazzitelli 
(1997) models. We assume an intrinsic cluster distance modulus 
of 8.10 and a reddening E(R - I) = 0.043. The dashed box in- 
dicates our estimate of the position and uncertainty of the LDB. 
The dotted contours show isochrones (derived from the Chabrier 
& Baraffe 1997 models) corresponding to the location of the LDB 
in this diagram. 



been depleted and we will use this to define the location of 
the LDB in the next section. 



5 THE LITHIUM DEPLETION BOUNDARY 
AGE OF NGC 2547 

5.1 The location of the LDB 

Figures E] and show the targets plotted in the I versus 
R — I, I versus I — Z and K„ versus I — K s colour-magnitude 
diagrams (CMDs). Objects classified as Li-rich members, RV 
members with no Li, non-members on the basis of their RV 
or Li-rich objects with discrepant RVs are represented with 
different symbols. There is considerable scatter in these dia- 
grams with no clear PMS locus visible - the isochronal ages 
discussed later in this paper are derived from stars of higher 
mass (see section 6.3.2). There are probably several sources 
of this scatter including binarity, rapid rotation, variability 
due to magnetic activity and photometric uncertainties. 

There is however a clear pattern among the Li abun- 
dances. The brighter objects are Li-poor - all RV cluster 
members with I < 17.22 have no Li. The fainter objects 
are Li-rich - all RV cluster members with / > 17.57 have 



Li. However, there are both Li-poor and Li-rich objects be- 
tween these two limits which potentially confuse the location 
of the LDB. Three of the four Li-rich objects in this region 
have broad Li lines suggesting rapid rotation which could re- 
sult in them appearing brighter and redder than suggested 
by a PMS isochrone for slowly rotating stars (Pinsonneault 
et al. 1998). Alternatively their location in all three CMDs 
suggests that these are unresolved binary systems, which 
constitute roughly 30 per cent of the cluster population at 
these colours (Jeffries et al. 2004). Perhaps a clearer way of 
looking at the data is to note that while there may be several 
reasons why a star with Li lies above the single star locus in 
the CMD, there are none (beyond photometric error) that 
would cause a star without Li to lie below the single star 
locus. The faintest object without Li has / = 17.49 and two 
others have I — 17.40. This suggests that the LDB must be 
at I ^ 17.4. We contend that the LDB of NGC 2547 is at 
/ = 17.54 ± 0.14 and at colours of R - I = 1.77 ± 0.03 and 
I — Z — 0.53 ± 0.02. Similar considerations lead us to place 
the LDB at K s = 14.86 ± 0.12 and I - K s = 2.70 ± 0.10. 
These locations are represented by the boxes in Figs. |S| 
and0 
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Figure 6. I, I-Z colour-magnitude diagram. Symbols are as in 
Fig. EI 



5.2 Age estimates 

The LDB location denned above can now be used to esti- 
mate the age of NGC 2547. The approach we take is the 
same as that described in Oliveira et al. (2003). Briefly, we 
use empirical relations between bolometric correction (for 7 
or K s ) and colour (R — I or I—K s ), along with the appropri- 
ate distance modulus, extinction and reddening (see section 
1 and below) to calculate the absolute I or K 3 magnitude of 
the LDB at any age, where we assume that the LDB repre- 
sents the point at which 99 per cent of the initial Li has been 
depleted and hence the EW of the Li I feature has dropped 
below about 0.2A. The age estimate can be performed for 
different evolutionary models. We note that our K 3 magni- 
tudes are on the 2MASS system and were converted to the 
CIT system (see Carpenter 2001) for the purposes of these 
calculations. 

Uncertainties in the derived age are produced by several 
factors (see Jeffries & Nay lor 2001). The ones we include in 
our analysis are: (1) The uncertainty in the position of the 
LDB in the CMD. These were estimated in the previous 
section and probably represent uncertainties of more than 1 
sigma. (2) Systematic errors in the photometric calibration. 
These are important for I, R — I and I — K 3 . Jeffries et al. 
(2004) estimated external uncertainties for stars with R — 
I > 1.5 of up to ±0.10 mag in I and R — I and hence 
there will be a systematic uncertainty of ±0.10 in the I — 
K s colours also. We assume a negligible external error for 
the K s magnitudes. (3) We assume an uncertainty in the 
reddening E(B - V) = 0.06 ± 0.02 (Claria 1982) and hence 
that E(R-I) = 0.043±0.014, E{I~K 3 ) = 0.092±0.031 and 
extinctions of Ai = 0.112 ± 0.037 and A K = 0.021 ± 0.007 
(Rieke & Lebofsky 1985). (4) We assume an uncertainty in 
the intrinsic distance modulus of 8.10±0.10 mag (see Jeffries 
& Tolley 1998; Naylor et al. 2002). 

Table [3] list the results of these estimates for the so- 
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Figure 7. K 3 ,I — K B colour-magnitude diagram. Symbols are as 
in Fig. E] The solid and dashed lines are 31.6 Myr and 20.0 Myr 
isochrones with colours and magnitudes taken directly from the 
models of Baraffe et al. (2002) and assuming an intrinsic distance 
modulus of 8.10 and E(I—K 3 ) = 0.092. The dotted contours show 
isochrones (derived from the Chabrier & Baraffe 1997 models) 
corresponding to the location of the LDB in this diagram. 



lar metallicity evolutionary models of Chabrier & Baraffe 
(1997), D'Antona & Mazzitelli (1997), Siess, Dufour & 
Forestini (2000) and Burke et al. (2004). The uncertain- 
ties quoted include all the effects discussed above, which 
are added in quadrature 2 . Table |3] also gives the bolometric 
magnitude of the LDB, and the masses and Tiff at the LDB 
for each model. 

As a comparison, we have also recalculated the LDB 
age of IC 2391 using our models and bolometric corrections. 
The LDB location was found by Barrado y Navascues et al. 
(2004) to be at I = 16.22 ± 0.10, R- I = 1.90 ± 0.05 and 
K 3 = 13.49±0.10, I-Ks = 2.75±0.10. An intrinsic distance 
modulus of 5.95 ± 0.10 and reddening equal to that of NGC 
2547 were also assumed. The LDB ages and parameters for 
IC 2391 are also given in Tabic [|] 

In addition to the random error bars quoted on the 
ages in Table |3] there are other small systematic uncertain- 
ties affecting the LDB age scale as a whole - the gravity 
dependence of the bolometric correction and the exact level 
of Li depletion that the LDB corresponds to (see Burke et 
al. 2004). The theoretical models of Baraffe et al. (2002) 
yield J-band bolometric corrections as a function of colour 
and age (and hence gravity). It seems that the relationship 
between bolometric correction and R — I is more age sen- 
sitive, with bolometric corrections lower by about 0.05 mag 



2 This assumes the sources of uncertainty are independent. This 
is unlikely to be true for the distance modulus and reddening 
which arc correlated in the sense that a larger assumed reddening 
leads to a larger distance modulus via main sequence fitting and 
both lead to a smaller LDB age. 
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for stars at 30 Myr compared with 5 Gyr. The effect is lim- 
ited to about 0.02 mag in the I — K s relation. Because the 
bolometric corrections we use are derived from high gravity 
field dwarfs a small adjustment to the bolometric corrections 
may be required. However, changes of this size would lead 
to LDB ages that were younger by only ~ f Myr. Similarly, 
Li depletion is so rapid that whether the LDB corresponds 
to 90 per cent depletion or 99.9 per cent depletion changes 
the LDB ages by only ±1 Myr. 

The ages in Table El show very close agreement both 
in terms of the ages using bolometric corrections deduced 
from the R — I and / — K s colours and between ages derived 
from different models incorporating different physics. Tak- 
ing the average result from both colours, the age of NGC 
2547 ranges from 34 Myr to 36 Myr and the age of IC 2391 
from 48 Myr to 53 Myr. The youngest ages in both cases are 
given by the D'Antona & Mazzitelli (1997) models. The dif- 
ference in LDB age between IC 2391 and NGC 2547 is very 
consistent at 14-17 Myr. 

As a final estimate of the LDB age we compared the 
absolute I and K s magnitudes at the LDB directly with 
the models of Baraffe et al. (2002), remembering to convert 
the 2MASS magnitudes to the CIT system predicted by the 
models. In other words we are using theoretical bolometric 
corrections produced by the cool model atmospheres. The 
results for NGC 2547 and IC 2391 are listed in the last col- 
umn of Table 01 and are very close to the estimates based on 
empirical bolometric corrections and the same evolutionary 
models (those of Chabrier & Baraffe 1997). 



6 DISCUSSION 

6.1 Comparison with Oliveira et al. (2003) 

The analysis of Oliveira et al. (2003) arrived at a slightly 
older lower limit to the LDB age for NGC 2547 than we 
have found in this paper. Oliveira et al. found no evidence 
for a population of Li-rich objects with / < 17.2, although 
they did find a few examples of Li detections of brighter 
stars which might be either binary members of the cluster 
or perhaps not members of the cluster at all - see below. At 
fainter magnitudes there were a few tentative Li detections, 
but to place limits on the location of the LDB, Oliveira et 
al. averaged the spectra of targets with 17.3 < I < 17.8 
and with I 17.8 separately. Evidence for an Li line was 
found in the average spectrum of the fainter sample, but 
not the brighter. Hence Oliveira et al. placed the LDB at 
17.8 < I < 18.3 and derived an age at least 3 Myr older 
than in this paper. 

With the benefit of hindsight and in the light of the 
much better data obtained with VLT/GIRAFFE we can ex- 
plain the probable cause of this discrepancy. We have found 
that the LDB is at / = 17.54 ± 0.14, but 5 of the 10 objects 
in the 17.3 < I < 17.8 sample of Oliveira et al. are brighter 
than this and hence unlikely to be Li rich. In addition we see 
from Fig. |S] that some level of contamination of the sample 
is also likely, such that perhaps 1-2 of the 5 objects with 
17.54 < I < 17.8 might be non-members and hence Li-poor. 
Given this, it is not hard to see why Li may not have been 
detected in the average spectrum of this subsample. On the 
other hand, the fainter sample should all be Li-rich , apart 



from contamination by Li-poor non-members at the level of 
^ 25 per cent, and hence an Li detection would be expected. 

There are also some discrepancies between the strengths 
of Li detected in this paper and those measured in Oliveira 
et al. (2003). We have observed three of the stars claimed 
to be Li-rich by Oliveira et al. Of these, we also detect Li in 
star 31 (= Oliveira no. 75, 075) and star 42 (= 091). But, 
we find a very low upper limit to the Li EW of star 50 (= 
O100), where Oliveira et al. claim a significant detection. 
We have also observed six stars with a significant Li de- 
tection - stars 5(=046), 27(=073), 39(=085), 41(=089), 
45(=92) and 53(=O106) for which Oliveira et al. found up- 
per limits to the Li EW. Of these, two have a detection 
consistent with the previous upper limit (stars 27 and 53), 
but the other four have Li EWs greater than the previous 2- 
sigma upper limit. Three of these (stars 39, 41 and 45) were 
classed as non-members by Oliveira et al. on the basis that 
their narrow band spectroscopic indices (the same as those 
used in this paper) were too low for stars of their colour. It 
now seems clear that the "non-members" found by Oliveira 
et al. were almost certainly instances of under-subtracted 
sky, because even the contaminants in the photometrically 
selected cluster candidates have spectroscopic indices sim- 
ilar to those of cluster members. An under-subtracted sky 
will result in lower spectroscopic indices, dilution of any Li 
EW and an erroneously optimistic assessment of the SNR 
of the spectrum. Hence it is not surprising that stars classed 
as non-members in Oliveira et al. do not have Li detections. 

6.2 Cluster membership and contamination 

The work of Jeffries et al. (2004) in determining the mass 
function of NGC 2547 relied solely on photometric criteria 
to select cluster candidates. The spectroscopic results pre- 
sented here offer an opportunity to check those photomet- 
ric criteria and to estimate the amount of contamination 
present in a photometrically selected sample. 

Figure [8] is an / versus R — I CMD showing the cata- 
logue of stars from which Jeffries et al. (2004) selected cluster 
candidates, their photometrically selected cluster candidates 
and symbols representing the status of stars that are anal- 
ysed in this paper. Note that the cluster candidates were 
also selected on the basis of their positions in the I versus 
I — Z and R — I versus I — Z diagrams (see Jeffries at al. 
for details) . There are several points arising from Fig. |H] 

(1) Looking at the distribution of RV selected members 
and RV non-members from this paper and comparing them 
with the distribution of photometrically selected candidates 
it appears that the photometric selection criteria employed 
by Jeffries et al. (2004) were quite appropriate. The presence 
of a "line" of five RV non-members just below the lower 
boundary of the cluster candidate distribution suggests that 
this boundary is well chosen, although we have also found 
two RV members (one Li-rich) near this boundary which 
were not photometrically selected as cluster candidates. 

(2) The upper boundary of the photometric selection is 
less well explored by our spectroscopic data. However this 
has lesser implications for the completeness and contamina- 
tion of the photometrically selected sample because the den- 
sity of stars in the CMD decreases rapidly towards brighter 
magnitudes. 

(3) Looking at spectroscopically observed stars which 



© 2004 RAS, MNRAS OOO.miTTI 



12 R.D. Jeffries and J.M. Oliveira 



Table 3. The LDB ages for NGC 2547 and IC 2391 using the / vs _R - / and K vs I - K CMDs using different models and a summary 
of model parameters at the LDB. 







Chabrier Sz 


D Antona & 


Sicss ct al. 


Burke ct al. 


Baraffe et al. 






Baraffe (1997) 


Mazzitelli (1997) 


(2000) z = 0.02 


(2004) 


(2002) theoretical 


NGC 2547 














LDB ages 


I vs R - I 


35.4 ±3.3 


34.0 ±3.7 


36.2 ± 3.4 


36.1 ±4.4 


36.3 ± 2.8 


(Myr) 


K s vsl - K s 


34.4 ±2.7 


33.8 ±2.8 


35.4 ±2.6 


35.4 ±3.1 


35.0 ± 2.6 


Parameters 


M bo i l 


9.58 ± 0.16 


9.58 ±0.16 


9.58 ± 0.16 


9.58 ±0.16 


9.63 ±0.16 


at the LDB 


Mass 1 (M©) 


0.17 ±0.02 


0.17 ±0.02 


0.15 ± 0.02 


0.17 ±0.01 


0.17 ±0.02 




Trff 1 (K) 


3250 ± 30 


3140 ± 20 


3230 ± 40 


3240 ± 40 


3240 ± 30 


IC 2391 














LDB ages 


I vs R - 1 


49.1 ±4.9 


47.3 ±4.3 


49.6 ±4.8 


52.5 ±5.7 


50.8 ± 3.9 


(Myr) 


K s vs I — K s 


50.4 ±3.8 


48.5 ±3.2 


50.5 ± 3.5 


54.1 ±4.4 


49.6 ± 3.8 


Parameters 


Mtot 1 


10.34 ±0.14 


10.34 ±0.14 


10.34 ±0.14 


10.34 ±0.14 


10.36 ±0.14 


at the LDB 


Mass 1 (M ) 


0.12 ± 0.01 


0.12 ±0.01 


0.11 ± 0.01 


0.12 ±0.01 


0.12 ±0.01 




Teff 1 (K) 


3140 ± 30 


3070 ± 20 


3100 ± 30 


3150 ± 40 


3130 ± 30 



1 - these are averages from the very similar results for the two diagrams 




Figure 8. I,R-I colour-magnitude diagram for our targets super- 
imposed on the general population towards NGC 2547, showing 
those objects identified as photometric cluster candidates by Jef- 
fries et al. (2004). The solid line shows the empirically calibrated 
30 Myr isochrones used to select members, which were generated 
from the models of Chabrier & Baraffe (1997 - solid line) and 
D'Antona & Mazzitelli (1997 - dashed line). Photometric selec- 
tion also made use of the I vs I — Z and R— I vs I — Z diagrams. 



are also photometrically selected candidates we find four 
RV non-members (plus another RV non-member which is 
Li- rich and a probable cluster binary - star 60) . In contrast 
we have found 50 RV members among the photometrically 
selected cluster candidates of which perhaps only 1-2 are 
random contaminants (none of which will be Li-rich). Hence 



the level of contamination among the photometrically se- 
lected cluster candidates is only about 10-15 per cent for 
the central regions of the cluster in the magnitude inter- 
val 16.5 < I < 19. Given that we have observed nearly all 
the cluster candidates within the VLT/GIRAFFE field of 
view, this translates into a spatial density of about 30 ± 15 
contaminants per square degree in this magnitude range, or 
45 ± 18 if we assume there are an additional 2 contaminants 
among the Li-poor RV members. This compares reasonably 
well, though is a little lower than the estimates from Jeffries 
et al. (2004) of 80 per square degree based on a population 
synthesis or the 64 ±26 per square degree based on the num- 
ber of photometric cluster candidates found in regions offset 
from the cluster. Differences in contamination levels of this 
size would have no significant effect on the mass and lumi- 
nosity functions for NGC 2547 determined by Jeffries et al. 
(2004). 

(4) There is scope to improve the precision of the LDB 
location in the CMDs. There are a further 8 photometric 
cluster candidates that we have not observed that lie inside 
the LDB error box in the I versus R — I CMD. However, 
the contributions to the error budget from uncertainties in 
the cluster distance and uncertain photometric calibrations 
for very cool stars play an equal role in the precision of the 
LDB age estimate. 

(5) Star 5 is Li-rich but quite clearly not a photomet- 
ric member of the cluster in any of the CMDs and not an 
RV member either. This star is too far above the cluster 
locus to be a a binary member. A single star of this colour 
would have to be younger than 50 Myr to have preserved its 
Li and such stars should be quite rare among the general 
field population. One possibility is that star 5 is a low-mass 
member of the Vela OB2 association which is sparsely spread 
over a ~ 10 degree diameter and which encompasses NGC 
2547. Vela OB2 is at a distance of ~ 400 pc (de Zeeuw et 
al. 1999) and has an age of a few Myr. A population of low- 
mass, Li-rich PMS stars has been found in this association 
by Pozzo et al. (2000, 2001). If star 5 were at a distance 
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of 400 pc, then its position in the CMDs suggests an age 
of about 4Myr (Baraffe et al. 2002), in excellent agreement 
with this hypothesis. Furthermore, star 5 has an RV which 
is 7.3 km s -1 larger than NGC 2547 average, which itself has 
a mean heliocentric RV of +12.8 ± l.Okms -1 (Jeffries et al. 
2000). Thus the heliocentric RV of star 5 is 20.1 ± l.Okms -1 , 
very close to the heliocentric RV of 18kms _1 obtained for 
the Vela OB2 PMS population by Pozzo et al. (2001). A 
similar explanation is likely for the Li-rich star RX2 found 
in the NGC 2547 study of Jeffries et al. (2003). It lies 2 mag 
above the NGC 2547 PMS locus, exhibits undepleted Li and 
has an RV consistent with either NGC 2547 or Vela OB2. It 
is fortunate for membership studies in NGC 2547 that the 
distance/age combination of Vela OB2 places its low-mass 
population well above the NGC 2547 locus in CMDs. 



6.3 Comparison of LDB ages with isochronal ages 

The LDB ages determined for NGC 2547 and IC 2391 in 
Table |S| (and also the Alpha Per and Pleiades clusters - see 
Jeffries & Naylor 2001) are very precise thanks to the steep 
luminosity dependence of the age at which Li is burned in 
fully convective stars. These LDB ages are also likely to be 
very accurate in an absolute sense. The range of models 
we have employed use a wide variety of equations of state, 
radiative opacities, atmospheres and convective treatments 
and yet the difference between the largest and smallest age 
estimate in the case of NGC 2547 is only 1.9 Myr! Of course 
there is always a concern that some physics neglected by all 
the models may be important, such as rotation or the pres- 
ence of supporting magnetic fields in the convection zone. 
Preliminary work suggests that the effects of these will be 
to increase the LDB ages, but only by a very small amount 
(Burke et al. 2004; F. D'Antona private communication). As 
the LDB ages are both precise and accurate they offer the 
opportunity to test the less well constrained physics that 
controls age determinations from higher mass stars. 



6. 3. 1 Nuclear turn off ages 

The LDB results for IC 2391, the Alpha Per and Pleiades 
clusters have shown that nuclear turn-off ages deduced from 
high mass stellar models without convective core overshoot- 
ing are younger by factors of 1.5 than LDB ages (Barrado y 
Navascues et al. 2004) . Unfortunately NGC 2547 may or may 
not be consistent with this picture. The nuclear turn-off age 
of NGC 2547 is ill-determined (55 ± 25 Myr based on mod- 
els with a small amount of convective overshoot - Jeffries & 
Tolley 1998), mainly because it hinges on the photometry 
and reddening of just one star (HD 68478, not HD 68468 as 
stated in Jeffries & Tolley 1998). 



6.3.2 Empirical low-mass isochrones in colour-magnitude 
diagrams 

The ages of young clusters can also be derived rather pre- 
cisely (but not necessarily accurately - see below) by fitting 
low-mass isochrones to cluster members in CMDs. The rate 
of descent towards the ZAMS is mass dependent but also 
depends somewhat on the treatment of convection and the 




Figure 9. An example V versus V — I CMD for presumed single, 
spectroscopically confirmed members of NGC 2547 (spots). Em- 
pirically calibrated isochrones at 20, 30 and 40 Myr arc shown, 
derived from the theoretical models of Baraffe et al. (2002, with 
a mixing length parameter of 1.9 pressure scale heights) with a 
distance, reddening and extinction appropriate for NGC 2547. 



Table 4. Empirical isochronal ages (in Myr) for NGC 2547 
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details of the stellar atmospheres. Here NGC 2547 is an ex- 
cellent cluster to check the consistencies of the age determi- 
nations. It has a well studied lower main sequence with very 
precise photometry and a well defined PMS locus. Unfortu- 
nately the absolute accuracy of the method is also depen- 
dent on how the model temperatures and luminosities are 
converted into colours and magnitudes before comparison 
with the data (or vice- versa). 

In previous papers (e.g. Naylor et al. 2002; Oliveira et 
al. 2003; Jeffries et al. 2004) a method was described for em- 
pirically deriving the relationship between colour and T e s, 
using the well-observed Pleiades cluster as a calibrator. We 
assume the Pleiades distance and age are 132 pc and 120 Myr 
(roughly the LDB age) and adopt the same relationships be- 
tween bolometric correction and colour as were used in deriv- 
ing the LDB ages. Then, assuming the same colour- T e ff rela- 
tion applies to clusters as young as NGC 2547 (which is not 
contradicted by theoretical atmosphere models) , isochrones 
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can be generated to find a cluster age. The method is insen- 
sitive to the age assumed for the Pleiades and changing the 
Pleiades distance simply changes the distance to the fitted 
cluster by a similar factor (and hence the LDB age - see the 
discussion in Oliveira et al. 2003). 

Figure shows examples using the D'Antona & Mazz- 
itelli (1997) and Baraffe et al. (2002, with a mixing length 
parameter of 1.9 pressure scale heights) models for an age 
of 30 Myr. Other examples of fits to NGC 2547 and IC 2391 
covering stars at higher masses and in other CMDs can be 
found in Naylor et al. (2002), Oliveira et al. (2003) and Jef- 
fries et al. (2004). We have uniformly reanalysed the V ver- 
sus V - I and I versus R- I CMDs for NGC 2547 to esti- 
mate low-mass isochronal ages using the same models used 
to estimate the LDB ages. We have taken spectroscopically 
confirmed stars from Jeffries et al. (2000, 2003) and excluded 
those stars which were considered to be binary systems in 
an attempt to define a clean single star isochrone. For stars 
with 1.0 < V - I < 2.5 and 0.5 < R — I < 1.5, spectro- 
scopically selected members of NGC 2547 lie on a very well 
defined PMS locus in the V versus V — I and / versus R— I 
CMDs (e.g. see Figs. l9l and 111(1 . Warmer stars have already 
reached the ZAMS and offer no age discrimination. Cooler 
stars appear to show a more scattered PMS (see Figs. |K| 
and II It . due to a combination of variability and photomet- 
ric errors, but may also suffer from systematic photometric 
calibration uncertainties. 

Model isochrones were generated using the Pleiades 
as an empirical calibrator of the colour- T B s relation and 
then matched to the NGC 2547 data using an intrinsic dis- 
tance modulus of 8.1, (which is constrained by higher mass 
stars that have already reached the ZAMS) reddenings of 
E{V - I) = 0.077 and E{R - I) = 0.043 and the corre- 
sponding extinction values. The results are given in Table 3] 
and an illustrative example of the isochrone matching shown 
in Fig. [5] The ages deduced from the two CMDs are very 
consistent, but vary from about 25 Myr for the D'Antona 
& Mazzitelli (1997) models (a value independently arrived 
at by Stauffer et al. [2003] using a similar technique), up to 
about 35 Myr for Baraffe et al. (2002) models with a mix- 
ing length set to 1.0 pressure scale heights. The photometry 
available for the IC 2391 cluster is more scattered, but sug- 
gests an isochronal age that is 5 to 10 Myr older using any 
of these models. 

The relative isochronal ages from different models are 
very precise because the alternative isochrones are almost 
parallel in the colour ranges considered. There are however 
uncertainties which shift the ages from all models approx- 
imately equally. The main sources of experimental error in 
this technique are the uncertainties in the adopted distance 
modulus and reddening of the cluster. Assuming the same 
values used to derive the LDB ages we find =p5Myr for a 
±0.1 mag change in distance modulus. This is a change in 
the same sense, but larger in magnitude than for the LDB 
age, where we find about =pl.6Myr for a ±0.1 mag change 
in the distance modulus of NGC 2547. A ±0.02 mag change 
in E(B — V) results in a change in the derived intrinsic dis- 
tance modulus (from the ZAMS stars) of about ~ ±0.1, but 
does not significantly change the isochronal age estimates 
from the lower mass stars since the reddening vector is close 
to parallel with the PMS isochrones in the CMDs. Statis- 
tical errors in deciding which isochrone fits best are only 



about ±2 Myr for NGC 2547 where, for 1.0 < V - I < 2.5 
and 0.5 < R — I < 1.5 the stars lie almost like "beads on a 
string" . The scatter is greater at redder colours in NGC 2547 
and also in IC 2391 at all colours, where an extra ±5 Myr 
uncertainty in the isochronal age is warranted. 

The main source of systematic error is the compatibility 
of the photometric calibrations of the Pleiades calibration 
data and the younger cluster data and the possibility that 
the colour- T e ff relation changes for stars with lower gravity. 
These factors are much less important for stars at 30 Myr 
with R — I< 1.5 and V — I < 2.5, which have gravities quite 
similar to Pleiades stars of the same colour and which are 
calibrated by many photometric standards of appropriate 
colour. We attach less weight to empirical isochrone compar- 
isons for cooler stars, where gravities are lower than in the 
Pleiades, where the photometric calibration may be uncer- 
tain by ±0.1 and where there appears to be more scatter in 
the photometry of genuine cluster members. However, even 
with these uncertainties it is encouraging that the empirical 
30 Myr isochrones shown in Fig. |5] are certainly compatible 
with the data for the cooler stars observed in this paper. 

Stauffer et al. (2003) have shown that the deduced 
isochronal age may also depend to a certain extent upon 
which colour indices are used. In particular they find that 
young, magnetically active stars suffer from a "blue excess" 
with respect to older stars. This would make ages deter- 
mined from the V versus B — V CMD problematic. However, 
they also show that this problem is much less serious using 
the V versus V — I CMD upon which our empirical isochronal 
ages for NGC 2547 largely rest. In any case, by using the 
Pleiades as an empirical calibrator of the colour- T e a rela- 
tionship we potentially nullify these problems because the 
cool stars of NGC 2547 have very similar magnetic activity 
levels to those in the Pleiades (Jeffries & Tolley 1998). We 
also find that isochronal ages determined from the V versus 
V — I and / versus R — I CMDs are entirely consistent. The 
same may not be true when it comes to comparing theoret- 
ical isochrones with the NGC 2547 data, because the best 
theoretical atmospheres currently contain no contribution 
from chromospheres, spots or plages (see section 6.3.4). 

The LDB age for IC 2391 is larger than its low-mass 
isochronal age when determined using the same models, al- 
though the discrepancy is perhaps not very significant at 
about 12 ± 8 Myr. It was partly this and partly the overes- 
timate of the LDB age for NGC 2547 that led Oliveira et 
al. (2003) to speculate that LDB ages were systematically 
higher than isochronal ages in low-mass stars, perhaps indi- 
cating the need for new physics in the low-mass evolution- 
ary models. Our new and precisely determined LDB age for 
NGC 2547 is only larger by a few Myr than the isochronal 
ages. The largest discrepancy occurs for the D'Antona & 
Mazzitelli (1997) models, and is only 10 ±6 Myr. Discrepan- 
cies of this size could reasonably be accounted for by a small 
decrease in the cluster distance modulus of about its esti- 
mated error bar or could be due to small systematic changes 
in the colour- T c g relationship with age. 

6.3.3 The empirical temperature scale 

The near consistency achieved between the LDB and 
isochronal ages is almost independent of the choice of model. 
This insensitivity to the chosen model arises from the intrin- 
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Figure 10. The relationship between intrinsic R — I and T e ff 
derived by requiring evolutionary models to match Pleiades data 
at an age of 120 Myr and distance of 136 pc. The relationships for 
the Chabrier & Baraffe (1997 - CB97, with overshooting of 1.9 
pressure scale heights), D'Antona & Mazzitelli (1997 - DaM97) 
and Siess et al. (2000 - SDFOO, solar metallicity model) are 
shown. Also shown are the semi-empirical relationships found by 
Houdashelt et al. (2000 - HoudOO) and Bessell (1979 - Bess79) and 
the empirical determinations of Kirkpatrick et al. (1993 - Kirk93) 
and Leggett (1996 - Leg96). 

sic similarity of the isochrones generated by the various mod- 
els at this age, but also by our calibration technique of fixing 
the age of the Pleiades and hence deriving a slightly different 
relationship between colour and T e s for each model. There 
are some differences in the isochrone shapes for cool young 
stars (see for instance Fig.|SJ, but unfortunately these occur 
at colours where uncertainties in the present photometric 
calibrations make it impossible to judge which models best 
describe the data. 

In principle we could compare the required intrinsic 
colour-Teff relationship with empirical and semi-empirical 
determinations of T a g to see which model correctly describes 
this relationship. Figure [TUl shows this comparison for R — I 
versus T c s, the plot for V — I versus T e g is qualitatively 
very similar. Given the uncertainties in the empirical T e g 
determinations - typically 200 K and the possible ±0.1 mag 
uncertainties in the photometric colours of very red stars we 
would be hard pressed to choose between the models in this 
way. 

Finally, we can compare the T e g at which the models 
predict the LDB (listed in Tabled with the T e ff correspond- 
ing to the colour of the LDB from Fig. ^| The T e ff at the 
LDB ranges from 3250 K, for the Chabrier & Baraffe (1997) 
models (which use a detailed model atmosphere) to 3140 K 
for the D'Antona & Mazzitelli (1997) models (which use a 
grey atmosphere). The derived intrinsic {R—I)-T E ff relations 
in Fig. HUI give very similar T e g values for all the models of 
3185-3220 K at (R - I) = 1.73. Given the possible 0.1 mag 
uncertainties in R — I, corresponding to T c ff uncertainties of 



0.5 1 1.5 2 2.5 

R-I 

Figure 11. An / versus R — I CMD for presumed single, spectro- 
scopically confirmed members of NGC 2547 (spots) and the radial 
velocity members identified in this paper (squares) . Isochrones at 
20, 30 and 40 Myr are shown and taken directly from the theoreti- 
cal models of Baraffe et al. (2002, with a mixing length parameter 
of 1.9 pressure scale heights) with a distance, reddening and ex- 
tinction appropriate for NGC 2547. 



~ 150 K, all the models are satisfactory in the sense that the 
isochronal T e g at the colour of the LDB and the theoretical 
Tefi at the LDB are consistent. 



6.3.4 Comparison with theoretical low-mass isochrones 

A purer test of the consistency between the LDB and low- 
mass isochronal ages is available using the models of Baraffe 
et al. (2002). These models, based on the evolutionary code 
of Chabrier & Baraffe (1997), predict both Li-depletion and 
the photospheric colours and magnitudes using complex 
model atmospheres. 

Figure [7] includes theoretical isochrones adjusted for an 
intrinsic distance modulus of 8.1 mag and reddening of E(I— 
K 3 ) = 0.092. These suggest an age of 25 ± 5 Myr for NGC 
2547, but there are additional errors due to the systematic 
photometric uncertainty of ±0.1 in the / magnitudes of these 
cool stars, ±0.031 in E(I — K s ) and ±0.1 in the distance 
modulus. Hence the age from this method is 25 ± 7 Myr, 
again in reasonable agreement with the LDB age found from 
the same stars and evolutionary models. The K B versus I — 
K s diagram for stars in IC 2391 around its LDB suggests an 
age of 40 ± 10 Myr, also in reasonable agreement with the 
LDB age (see Fig. 9b of Barrado y Navascues et al. 2004) . 

The cool stars we have considered here have rather im- 
precise 2MASS near IR magnitudes. To improve on this we 
use the spectroscopically confirmed sample of brighter NGC 
2547 stars discussed in section 6.3.2. The comparisons be- 
tween theory and data are shown in Figs.[TTllIland[T31 Note 
that the Baraffe et al. (2002) colours and magnitudes on the 
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Figure 12. An / versus I-J (2MASS system) CMD for NGC 
2547. Symbols and lines are as for Fig. 1111 



I— J < 3. In all the CMDs, but particularly for I versus R—I 
there is a hint that the stars are not arrayed parallel to the 
isochrones even in these colour ranges. At redder colours this 
becomes more than a hint. The data and models are clearly 
discrepant at any age for R — I > 1.2 (a deficiency in the 
models already noted by Baraffe et al. [1998] and probably 
due to missing atmospheric molecular opacity in the opti- 
cal spectrum of cool objects) and a drift of stars towards 
younger ages for I — J > 1.5 is also apparent. However, re- 
call that systematic errors of order 0.1 mag are possible in 
the colours of the red stars so overall the theoretical NIR 
colours and models do a reasonably good job of modelling 
the NGC 2547 data at or just below the LDB age derived 
from the same model. It is interesting to note though that 
the theoretical model atmospheres do not simulate spots, 
plages, chromospheres or other manifestations of magnetic 
activity in these young stars. Stauffer et al. (2003) found 
that cool Pleiades stars seemed to show a small if-band ex- 
cess with respect to older stars. There is a hint of this also 
in our data, where the K versus I — K CMD gives a lower 
isochronal age than the other CMDs. This would be brought 
back into line by the addition of a ~ 0.1 mag I — K excess 
in the model atmospheres. 
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Figure 13. A K 3 versus I-K s (2MASS system) CMD for NGC 
2547. Symbols and lines are as for Fig. 1111 

CIT system have been converted to the 2MASS system for 
these plots. 

The agreement between the LDB age of NGC 2547 from 
the Chabrier & Baraffe (1997) models (35 ±3 Myr) and the 
theoretical isochrone fits from the same models is also very 
good. The theoretical I versus R—I isochrones, like the em- 
pirical isochrone modelling in Jeffries et al. (2004), suggest 
an age of 30 ± 5 Myr for R — I < 1.2, the / versus I — J 
plot suggests an age of (30 ± 5) Myr for I — J < 1.5, whilst 
the / versus I — K s CMD suggests an age of (25 ± 5) for 



7 CONCLUDING REMARKS 

Deep, intermediate resolution spectroscopy obtained with 
the VLT/GIRAFFE spectrograph has been used to study a 
group of 63 low-mass (0.08 — 0.3 Mq) candidate members of 
the young open cluster, NGC 2547. Radial velocity measure- 
ments and the presence of strong Li I 6708A absorption in 
the fainter stars, have confirmed membership for at least 50 
objects. We have hence validated previous work that selected 
low-mass members and estimated the extent of sample con- 
tamination based solely on photometric selection (Jeffries et 
al. 2004). The conclusions of that work regarding the mass 
function and luminosity function of NGC 2547 should be 
sound, at least down to the substellar boundary. 

From the sample of confirmed members we have deter- 
mined the magnitude at which lithium switches from being 
completely depleted, to being found close to its initial abun- 
dance, as 7 = 17.54 ±0.14, K s = 14.86 ±0.12 and at colours 
of R- I = 1.77 ± 0.03 and I — K s — 2.70 ± 0.12. 

Using several PMS evolutionary models which incorpo- 
rate a variety of physical approximations we find an almost 
model-independent lithium depletion boundary age (LDB) 
of 34-36 Myr, with a precision of 10 per cent. This LDB age 
is only slightly larger than the ages of 25-35 (±5) Myr de- 
duced from isochrone fits to stars from O.O8-1.2M0 using 
the same models and empirically calibrated relationships 
between colour, bolometric correction and T e a. It is also 
slightly larger than ages of 25-30 (±5) Myr deduced from 
isochrone fits in the 1,1 — J and 1,1 — K 3 CMDs using 
only the theoretical colours and magnitudes of Baraffe et al. 
(2002). However, perfect agreement with the LDB ages can 
be obtained with a distance reduction to the cluster of only 
0.1 mag, approximately equal to its estimated error bar. 

On the basis that the isochrone fits to low-mass are com- 
patible with LDB ages derived from the same models we find 
that all of the PMS models we have considered (Chabrier 
& Baraffe 1997; D'Antona & Mazzitelli 1997; Siess et al. 
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2000, Baraffe et al. 2002) give acceptable results. This is de- 
spite possible systematic uncertainties regarding the intrin- 
sic colours of cool, magnetically active young stars (Stauffer 
et al. 2003). At the age of NGC 2547, significant differences 
in the isochrones predicted by these models occur at colours 
so red that uncertainties in photometric calibration prevent 
any decisive test between them. However, this result now 
seems to remove any pressing need to introduce new physics 
such as rotation or extra magnetic support into the PMS 
models and lends confidence to isochronal ages determined 
(using the methods outlined here) from low-mass stars be- 
tween 30 and 120 Myr. This is an important result, because 
PMS isochronal ages can easily be determined from photom- 
etry obtained on small telescopes, whereas even for nearby 
clusters, determining LDB ages stretches the capabilities of 
8-m telescopes. 

The differences in the isochrones predicted by the vari- 
ous flavours of model, arising from variations in their treat- 
ment of convection, atmospheres, equations of state and in- 
terior opacities, become more apparent at younger ages. This 
work demonstrates that finding the LDB is certainly possi- 
ble down to I ~ 18.5 and therefore the technique could and 
should be applied to more distant clusters in the 10-20 Myr 
age range. 
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